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1 Introduction 



The Kadomtsev-Petviashvili (KP) equation and the whole KP hierarchy as 
well as the multicomponent KP and two-dimensional Toda lattice hierar- 
chies are among the most important ingredients of the theory of integrable 
partial differential equations. They arise in various fields of physics from 
hydrodynamics to string theory. 

Integrable hierarchies have been described and studied within the frame- 
work of different approaches. The Sato approach and the 5-dressing 
method (see 0- and review |jlO|) are among the most beautiful and pow- 
erful among them. These two approaches look like completely different. The 
Sato approach is in essence algebraic with its infinite-dimensional Grassman- 
nian, pseudo-differential operators, Hirota bilinear identity and r-function. 
In contrast, the 9-dressing method, based on the nonlocal 9-problem, mostly 
uses the technique of the complex variable theory. 

Recently we proposed {^^, [0) an approach which combines the char- 
acteristic features of both these methods, namely the Hirota bilinear iden- 
tity from the Sato approach and analytic properties of the wavefunctions 
from the 9-dressing method. This analytic-bilinear approach is based on the 
generalized Hirota bilinear identity for the wave function ilj{X,fi) with sim- 
ple analytic properties (Cauchy-Baker-Akhiezer (CBA) function). The CBA 
function (or, to be more precise, its 'nonsingular part' x('^)/^) is defined in 
some region G of the complex plane and it is characterized by the generalized 
Hirota bilinear identity on the boundary of this region [l^. The CBA func- 
tion arises both in 9-dressing method ^ and in finite gap integration method 
[Q, and it can be introduced also to the context of general Grassmannian 
approach 0], 

According to generalized Hirota bilinear identity, the CBA function is a 
functional of the function g defined on the boundary of G (in fact, g is a map 
of dG to (i.e. C-0), or, in other words, it is an element of loop group). One 
can choose some special subgroups of the generic loop group parametrized by 
finite number of variables ('times') to get particular examples of integrable 
equations. This is a standard approach in the 9-dressing method. One may 
also try to use some generic loop group and extract a full infinite set of 
equations. So the notion of the hierarchy arises. The dependence of the 
CBA function on g encodes the dynamics of the system. The hierarchy can 
be written down in variational form, i.e. as equations describing variations 
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of the objects of theory under the arbitrary variation of group element g. To 
get an hierarchy in terms of times, one should introduce a parametrization 
of the loop group. 

The distinctive feature of analytic-bilinear approach is that our primary 
object is the CBA function and hierarchies are written in the form of equa- 
tions for the CBA function (we call these equations generalized hierarchies). 
The analytic-bilinear approach allows us to describe generalized hierarchies in 
a compact finite form. Such a generalized hierarchy contains basic integrable 
hierarchy, hierarchy of equation for wave functions (modified hierarchy) and 
hierarchy of equations for CBA-type wave functions (we call it singularity 
manifold equations (SM) hierarchy, because in the scalar KP case it coincides 
with the singularity manifold equations obtained from the Painleve analy- 
sis of KP; however, our contexts is completely different), and corresponding 
linear problems. The relation between three hierarchies of equations can be 



described in terms of invariants of Combescure symmetry group ( |Tl 

The generalized Hirota bilinear identity provides various functional equa- 
tions for the CBA function. The resolution of some of them leads to the 
introduction of the r-function, while others give rise to the addition theo- 
rems for the r-function. The properties of the r-function, such as the associ- 
ated closed 1-form and its global definition are also arise in a simple manner 
within the analytic-bilinear approach. 



The first part of this paper was devoted to the scalar KP hierarchy. 



In this case G is just a domain of the complex plane. Without the loss of 
generality, we choose G as unit disc. The group element g is defined on the 
unit circle and can be represented in a form 

^(A) = A"^in X ^out 

where n is the index of the function, gi^ and gout are functions analytic and 
having no zeroes respectively inside (outside) the unit circle on the complex 
plane (i.e. elements of the loop groups r_ and r_|_). In the scalar KP hi- 
erarchy case, we have no opportunity to treat index variable without some 
additional construction, so we put index to zero. The functions ^^out and gin 
are parametrized in a standard way by infinite sets of 'positive' and 'negative' 
times. But it appears that in the scalar KP hierarchy case dependence of the 
objects of theory on g^^i (i.e. on negative times in our notations) is trivial, so 
the scalar KP hierarchy is defined effectively by the subgroup r+ of the loop 
group F. 
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In the present second part of the paper we consider the generahzed multi- 
component KP hierarchy and the two-dimensional Toda lattice hierarchy (for 
other approaches to the usual multicomponent KP hierarchy and 2D Toda 
lattice hierarchy see i, 0, §, §, p|, [|16l). 

There are two equivalent approaches to investigate mult i- component hi- 
erarchy. One could start with the generalized Hirota identity for the matrix- 
valued CBA function. In this case G is again a unit disc, and 5^ is a diagonal 
matrix function on the unit circle (i.e. a set of N functions Qn)- Another way 
is to start with the scalar CBA function defined in the region G consisting of 
N disconnected domains (unit discs for simplicity) . In this case the dynamics 
is also defined by the set of functions and the CBA function consists, 
in fact, oi N X N disconnected components. These two representations are 
exactly equivalent and are associated with the multicomponent loop group 
r^. In our paper we will use both of them. 

As in the scalar case, the objects of multicomponent KP theory depend 
trivially on the subgroup T^. The difference with the scalar case consists in 
the presence of discrete variables (indices), though only the lattice vectors 
with zero general index are admissible. 

Since the generalized Hirota bilinear identity is our primary object, it is 
quite natural to proceed from the hierarchy corresponding to the simplest 
choice of G - the unit disk - to the hierarchies corresponding to G consisting 
of several unit discs, i.e. from the scalar KP hierarchy to multicomponent 
KP hierarchy. But on the other hand, the loop group is also a fundamental 
object of the theory, so one would like to find an hierarchy with a full-scale 
action of generic loop group on the objects of theory. This hierarchy is 
known as Toda-lattice hierarchy. In our approach it can be introduced as a 
two-component KP hierarchy with the special choice of G consisting of two 
discs with the common boundary (the unit circle): one disc around origin and 
one disc around infinity. In this case, up to some trivial dynamics, one can 
glue two special loops defining the two-component KP hierarchy to obtain a 
generic loop. In the case of Toda lattice hierarchy the discrete index variable, 
positive and negative time variables get their natural interpretation. 

We first consider the multicomponent KP hierarchy. We start with the 
generalized Hirota identity for the matrix-valued CBA function and derive 
various finite functional equations for this function. Part of these functional 
equations present the compact form of the generalized multicomponent KP 
hierarchy which includes equations for potentials (the standard multicompo- 
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nent KP hierarchy), equations for the wave functions (the modified muhicom- 
ponent KP hierarchy), equations for CBA-type wave functions (singularity 
manifold equation hierarchy) and corresponding linear problems. 

Resolution of functional equations introduces the set of r-functions. We 
discuss also the associated closed 1-form and its global definition. The sim- 
plest integrable systems are given by the resonantly interacting wave equa- 
tions and by the scalar Darboux system which is its modified partner. 

We present a compact form of equations of multicomponent KP hierarchy, 
modified KP hierarchy and singularity manifold equation hierarchy, using 
shift operators. Similar formulae for the first two cases can be found in [|1^ . 

The two-component case which corresponds to the generalized Davey- 
Stewartson (DS) hierarchy is considered in detail. 

The description of the generalized Toda lattice hierarchy within the analytic- 
bilinear approach requires minor modifications of the multicomponent KP 
case. For simplicity we consider here only the Abelian Toda lattice hierar- 
chy. It appears that modified Toda lattice equation (i.e. equation for the 
Toda lattice wave functions) coincides with the 2D Volterra system, which 



was recently discussed in literature in different contexts |]18| , ||19[ . We derive 
2D Toda lattice singularity manifold equation, which was also discussed in 
The invariants of Combescure transformations group (|ll|, |T^) are 



found. We also present the hierarchies of equations using shift operators. 

The paper is organized as follows. The generalized Hirota bilinear identity 
and some basic formulae are discussed in section 2. Generalized multicompo- 
nent KP hierarchy is considered in section 3. Compact forms of equations of 
generalized multicomponent KP hierarchy are presented in section 4. Section 
5 is devoted to generalized DS hierarchy. The r-function for the DS hierar- 
chy is introduced in section 6. The r-function and corresponding addition 
formulae for the multicomponent KP hierarchy are considered in section 7. 
Closed variational 1-form and global definition for the r-function are given 
in section 8. Basic equations for the generalized 2D Toda lattice hierarchy 
are derived in section 9. Concrete equations associated with the generalized 
2DTL hierarchy are presented in section 10. Addition formulae for the 2DTL 
hierarchy are given in section 11. 
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2 Generalized Hirota identity, functional equa- 
tions and generalized multicomponent KP 
hierarchy 

The analytic-bilinear approach starts with the generalized Hirota bilinear 
identity 

Xii^, yu; 9i)9iii^)92^i^)xi>^, 92)dv = 0, (1) 

where xi^^ /^j 9) is the CBA function normalized by the condition x ~^ ^ 
as X —>■ fi and /x) is analytic function of both variables A, for A 7^ /x 
in the domain G with the boundary dG. The quantity g = g{X, x) specifies 
the dependence on the infinite set of dynamical variables x = (xi , X2 , ...)• Iii 
terms of the wavefunction ip{X, fi; g) defined via 

V'lA,/^;^) = ^~^(/i)x(A,/i;^)^(A) (2) 

the generalized Hirota bilinear identity looks like 

/ '^{^,fJ';9i)'^{>^,^;92)du = 0. (3) 

JdG 

Note that the generalized Hirota bilinear identity is bilocal both in dynamical 
and spectral variables. 

In this paper we consider the case for which x and ip take their value in 
the full linear group GL(N,C), G is a unit disc with a center A = 0. Elements 
g{X,x) belong to maximal abelian subgroup of GL(N,C). 

It is important to note that there is an alternative way to treat the matrix 
problem (|l|), starting from the sca/ar identity but with the region G consisting 
of several disconnected domains. In this case one has to choose domain G 
consisting of N unit discs with the center at A = 0. Now the scalar CBA 
function effectively consists of N x N disconnected components and can be 
interpreted as matrix function defined on the unit disc: 

Xa/3(A,/i) = Xe Dp, ne Da. (4) 

The function g in this framework is a set of functions defined on the unit 
circle. Thus, rewriting the scalar identity in matrix notations, we get exactly 
the same matrix generalized Hirota identity, where g is a diagonal matrix. 
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So the matrix structure in the scalar case is hidden in the multicompo- 
nent domain of definition of the CBA function. The advantage of the scalar 
treatment is the opportunity to preserve most of the structure of the one- 
component KP hierarchy case almost intact and to use geometric intuition in 
the matrix case. We will use both representations of the generalized Hirota 
identity in the paper. 

The dependence on dynamical variables is hidden in the set of functions 
g defined on the unit circle. The most natural choice would be to treat sets 
of arbitrary functions on the circle having no zeroes, forming the group by 
multiplication. Arbitrary g can be represented in the form 

where discrete variable n is the index of the function, g^^ and (^out are functions 
analytic and having no zeroes respectively inside (outside) the unit circle on 
the complex plane. We will use the following parametrization for g^y^ and 

gaont = exp X(„)„A"")j , (5) 

gain = exp 1^ ^ a;(„)„A"")j . (6) 

In these formulae we have omitted 'the zero time' Xq, that defines simply 
some constant factor in front of g and, according to (0), it doesn't change 
the function x- 

So, in the general case the set of functions g^ can be parametrized by 
A^ discrete times na and A^ infinite sets of continuous variables Xan, n & Z. 
But in the case of KP hierarchy the dependence of the function x(A, fi; g) on 
gia (in other words, on the 'negative times') is trivial; according to ([|) (the 
'scalar' form), we get 

X(A, fi;gx g^^) = gr^^{fj,)x{X, /x; ^)^m(A). (7) 

It occurs that other important objects (KP solution, r-function) also depend 
on negative times trivially. Thus, for simplicity, from the very beginning we 
omit the dependence on the negative times. 

As for discrete times, we will see, that the case of multicomponent KP 
hierarchy allows only sets of discrete times with the zero general index 
(Eliri. = 0). 



7 



So the case of multicomponent KP hierarchy seems a bit asymmetric with 
respect to positive and negative times, and it also doesn't allow to consider 
arbitrary lattice vectors of discrete times. The most natural way to resolve 
these problems is to consider multicomponent Toda-lattice hierarchy (see 
below). The objects of Toda-lattice hierarchy have a natural dependence on 
positive, negative and discrete times. 

To find special solutions for (P, one may use some additional construction 
(e.g. 9-dressing method or algebro-geometric technique). But it is possible 
to develop a procedure of solving (|l|) in terms of ([l|) itself. In this case 
the function x(A,/i,(y' = 1) = xo(A,/i) should be treated as initial data, 
and one can start from arbitrary xo{^, f^) possessing the necessary analytical 
properties. Then the dynamics of the function x(A, fi) can be defined via the 
variational form of identity (|l|) 

^x(A, /i; 5-) = f x{^, yu; 5^)— rYX(A, z^; g)du . (8) 

So the problem of solving ([^) is reduced to (infinite number in terms of 
variables x„) differential equations with initial condition. Taking the second 
variation in 6{lng), it is easy to prove that this system is self-consistent 
(i.e. the cross-derivatives are equal), and it is a correctly defined variational 
equation for the function x(A,/i, (?). Going further, one could take higher 
variational derivatives (all of them have a very simple form) and represent a 
formal solution in terms of Taylor-type series. 

Another approach to solving (|l]) is to extract from it some linear integral 
equations. Indeed, let us choose some initial condition xo(A,/i). Then for 
the function x(A, fi, g) from (|l|) we have a pair of equations 

x(i^,At;5')fi'('^)Xo(A, z/)rfz/ = 0, (9) 

dG 

Xoii^, fJ')9~^ii^)xiK i^; 9)diy = 0- (10) 

dG 

These relations can be interpreted as linear integral equations for the function 
x(A, fj,] g); integral equations of such type arise in the context of nonlocal Rie- 
mann problem. Using the technique developed in |13|, one may show that if 



both of these equations are solvable (in the class of functions xi^^^ /^) with the 
analytical properties specified above, here they play a role of normalization) , 
then the solution for both equations is the same and it is unique. 
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There is a very important class of functions g corresponding to some 
special discrete transformations - meromorphyc functions - for which the 
dynamics can be found explicitly. Namely,the solution for both equations 
(P), (p!OD exists (and it is unique), if and only if meromorphic function g has 
equal number of zeroes and poles in G (taking into account the multiplicity), 
and it is given by the determinant formula [0. This statement means that 
the dependence of x on discrete (index) variables is correctly defined for the 
zero general index case and, moreover, it can be represented in explicit form 

a- 



3 Generalized multicomponent KP hierarchy 

Proceeding with the multicomponent KP case, we choose the following parametriza- 
tion for g^nt (matrix representation) 



5((A,x) = exp 



N oo 



a=ln=l ^ 



(11) 



where the matrices Ea form a basis of the commutative subalgebra of diagonal 
matrices: {Ea)i3-y = SapSp^ (a, /5, 7 = 1, A^). So we have infinite series 
of dynamical variables X(^a)n- 

Equations (0) or (|^) with g of the form (|TT]) generically define the multi- 
component KP hierarchy. To extract from them finite functional equations 
similar to the scalar case we consider the Hirota identity (0) for the dynamical 
variables x'^^-) ^ and nrelated in a special way, namely 

/ 1 1 

^{a)n ~ ^ia)n = —a-a ~ ~K ^ (12) 

where aa, (0 > a > A^) are arbitrary parameters. We will also use vector 
notations x = (xi, xn), x^ = {x(^a)i, X(a)2, •••)• this terms the preceding 
formula looks like 

X(«) - X(„) = [aa] - [ba] , (13) 

or 

x'-x=[a]-[6], (14) 
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where [si]{a)i = jo^- For such x and x' due to the formula log(l— e) = J2'^=i ^ 
one has 

u - an 



u — h, 



- 0>a,p>N. 



(15) 



Substitution of the expression (15) into (|T|) gives 



— 



Xa/3(A,Ai,x+ [a]) 



A — a^' 
X-b 



Xa/3(A,/i,X+ [b]) + 



AT 



^(6^ - a^)x7/3(A,6T,,x+ [a])xa7(&7,/i,x+ [b]) = 0. (16) 

7=1 

Equation (ITBp represents the final functional form of the generalized multi- 
component KP hierarchy. It encodes the usual multicomponent KP hierarchy 
and its modified version together with their linear problems. To demonstrate 
this fact we, similar to the scalar case |T^, consider the particular form of 

= 0, > a > 

1 - — 1 Xa/3(A,/U,X+ [a]) - ( 1 - ^ ) Xa/3(A,/i,x) = 



equation (|l^) with 6q, = 0, 0>q;>A^ 



N 



«7X7/3(A, 0, X + [a])xa7(0, /i, x) = 0, 

7=1 

or, in terms of the CBA and Baker- Akhiezer functions 



(17) 



N 

ipa^iX, /i, X + [a]) - tpa^iX, = Y^ a^il)^p{X, x + [a])V^„^(/i, x)(18) 

7=1 

where V7/3(A,x) = x(A, 0, x)5(/3(A), V'7/3(/i,x) = x(0, yU, x)5(-^(/i). The expan- 
sion of ([18|) in Taylor series in a(^) gives rise in the first order to the following 
equation 

^V'a/3(A,/i,x) , (X ^-r ( ^ 

^ = ^"7/3 ( A, X)^a^ (/i, X) . (19) 

OX{^) 1 



Evaluating the equation ([T7| ) with 0:7^/37^77^0; at A = /i = 0, one obtains 
the system 



p P 

1 a-yi ■^p , 



(20) 
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where Pq,/? = Xa/3(0,0,x), = which is nothing but the Darboux 



system in terms of rotation coefficients [21| or resonantly interacting waves 



equation Equation ( jn\ ) with a 7^ 7 and /z = gives the system 

^ = P.,v^,, (21) 

where 

TV 

vl/„ = 5] irfA^(A)^„^(A,x) (22) 

(3=1 •' 

and //3(A) are arbitrary functions. The system (^TD is a well-known linear 
problem for equations (^). 

Equation (|l^ with /3 7^ 7 and A = gives the adjoint linear problem for 



the system 

^"^^ - P^,^, , (23) 



9^7 

where \E'/3 = I]a=i § djjfaifxjipaisil^, x) and /a(/i) are arbitrary functions. 
Further (^) implies that Pa-y = Substituting this expression into 

), one gets 

dpd,^^ = ^p~\d,^p){dp^^) + ^,~\dp^,){d,^^), (24) 

where djs = This is a well-known Darboux-Zakharov-Manakov system 
which has been discovered in the theory of systems of conjugate surfaces ffist 
[21] and then has been rediscovered within the 9-dressing method in [|^. This 
system is a modified partner of the system (|20|) . To get the linear problem 
for it, it is sufficient to substitute the expression Pay = -^^^ into (p^). 

Thus the system (|1^ contains the equation (pO]), its modified partner 
(E¥) and their linear problems. 

The higher order terms in the expansion of ( [T7| ) in give rise to the 



higher analogues of equation (p^Ql) . It is a straightforward check that these 
equations contain the higher order equations (|20|), their modified analogues 
and corresponding linear problems. 

The generalized multicomponent KP hierarchy described by functional 
equation (^) has rather special structure. In particular, it contains N in- 
dependent generalized scalar KP hierarchies. They are obviously given by 
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equations (|I6]) with a = P and = for all 7 except 7 = 0. It is also quite 
clear that the number of effective independent dynamical variables is equal 
to three. For instance, the systems PO] ) and ( P^ are overdetermined and 
represent a collection of commuting closed systems defined for every triad of 
dynamical variables. 

Note also that one can get the multicomponent KP hierarchy starting 
with the scalar generalized Hirota identity (0). In this case one has to choose 
disconnected domain G consisting of N unit discs Da with the center at 
A = 0. In this case the scalar CBA function effectively consists of N x N 
disconnected components and can be interpreted as matrix function. So the 
matrix structure in the scalar case is hidden in the multicomponent domain 
of definition of the CBA function. The function g{X) is defined on the set of 
unit circles Sa and consists in fact of A^ components ga defined on a unit 
circle. 

The advantage of the scalar treatment is the opportunity to preserve most 
of the structure of the one-component KP hierarchy case almost without 
modification. 

Using the parametrization @, we represent the function g{X) (defined on 
the set of A^ unit circles Sa) as 

TV / 00 \ 

^7(x,A)= nexp ' (25) 

a=l \i=l / 

where 

Ka{\) = j, AG D„; 

KaiX) = 0, \i D„. 

In this case the transformation group element g^ providing just the same 
shift of the x variables takes the form 

g,{v) = {g,{v)g,\v)) = ''-^ z/GD„. (26) 

U Oa 

The analogue of equation (|16]) reads 

9oifJ')xi\ /U, X + [a]) - go{\)xi\ yU, x + [b]) + 

N 

- a7)x(A,&(7),x+ [a])x(&7,/i,x+ [b]) = 0. (27) 

7=1 
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It seems to be a scalar equation, but if we take into account that the func- 
tion X consists of a set of disconnected components and denote Xa/si^, f^) = 
x(A, /i), A G Df3, fi G Da, then this equation gives rise exactly to the system 

The analogue of equation (|T7|) is respectively 

N N 

Y[il- aaKM)xiK X + [a]) - JJil- a„i^'c.(A))x(A, /i, x) = 

0=1 a=l 
N 

^a(^)X7(A,x+ [a])x7(/i,x) = 0, (28) 

7=1 

and in terms of CBA and BA functions one gets 

N 

ip{\ /i, X + [a]) - ip{\ /i, x) = ^ a^il)^{\, X + [a])V^^(/i, x), (29) 

7=1 

where we have introduced notations 

X7(A,x) = x(A, 0^,x), 

X7(/i,x) = x(0^,/i,x), (30) 
V^„(A,x) = 5((A)xa(A,x), 

V'a(/i,x) = ^~^(/i)Xa(/X,x), 

Oq, denotes the center point of the disc D^. 

We will also use following notations, most of them exactly corresponding 
to notations introduced in frame of matrix representation: 
for potentials 

Xa/3(x) = Xr(0;3,0a,x), 

Xr(A, yu, x) = x(A, /i, x) - 5ap—^, X e Da, fi e Dp 

A — 

for wave functions 

<l>(x)= // rfArf/i/(A)/(/i)^(A,/x,x), 
^,(x) = / dA/(A)^,(A,x), 

JdG 

^„(x) = / dfxf{fx)^IJail^,x.), (31) 

JdG 

where /(A),/(/i) are arbitrary functions defined on the boundary of G (i.e. 
on N unit circles). 
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4 Equations of generalized KP hierarchy in 
compact form 

The basic equation of the multicomponent KP hierarchy ([T7| ) looks very sim- 
ple and contains all the information about generalized hierarchy - hierarchy 
of linear problems and hierarchy of equations. However, sometimes this in- 
formation could be excessive and one would prefer to have a compact form 
of the equations of the hierarchy. As we have demonstrated in the previous 
section in the case of DS hierarchy, it is possible to extract all the hierarchy 
of equations from (|T^ in terms of shift operators. To derive the equations 
of the KP hierarchy, we will use generic shifts of times. For each triad of 
different shifts (i.e. shifts with different parameters) it is possible to derive 
a closed system of equations, representing, in fact, the multicomponent KP 
hierarchy. 

As we will see later, the multicomponent KP hierarchy in a general case in- 
corporates also discrete times not reducible to the shifts of continuous times. 
So equations of multicomponent KP hierarchy can also contain these es- 
sentially discrete times. It is possible to include them to equations of this 
section, but now we will derive the equations of KP hierarchy containing 
only continuous times. Discrete times will be discussed later in the context 
of Toda lattice hierarchy. 

Thus, we will derive a compact form of the equation of multicomponent 
KP hierarchy, modified KP hierarchy and singularity manifold equation hi- 
erarchy. It is interesting to note that these equations are obtained in fact in 
the form of lattice equations, and all the hierarchy is given by the expansion 
of these lattice equations with respect to the parameters of the lattice. 

In this section we will use the matrix form of the equation (0), i.e. 

Ax(A, /i) + \x{K f^)A - -A ■ Tx(A, fi) = x(0, fi)A Tx(A, 0). (32) 
A /i 

Here the common argument x is omitted, the function x(A,;U, x) is N x N 
matrix-valued function, A is a diagonal matrix, Aij = Sijai, and we have 
introduced shift operator and difference operator 

Txi>^, yLt, x) = x(A, /i, x + [a]), 
Ax(A, /i, x) = (T - l)x(A, /i, x). 
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If one takes the first term of expansion of the relation (^) at A = 0, /i = 0, 
one obtains an identity containing the potential u = Xr(0,0,x) and two 
higher terms of expansion of Xr(-^; x). It is possible to cancel these terms 
by combining three equations (|32D with different shift parameters a^, k G 
{1, 2, 3} and get an equation for the potential u. 

First, it is easy to cancel the terms containing Ai/X combining pairs of 
equations 

Aix{X, fJ^)Aj - Ajx(A, fi)Ai + -{AiTixiX, fi)Aj - AjTjxiX, fi)Ai) = 
X(0, fi)A,Ta{X, 0)A, - x(0, fi)A,T,x{X, 0)Ai . (33) 

And then, multiplying this equation by the operator AkTk from the left, 
taking sum by different permutations of indices i,j,k and using the first 
term of expansion by A, /i at A = 0, yU = 0, one obtains the equation 

eijkAkTk{AiU - uAiTiu)Aj = 0, (34) 

(ijk) 

where u = Xr(0, 0, x), summation goes over different permutations of indices. 
This extremely simple equation generates all equations of multicomponent 
KP hierarchy (not containing essentially discrete variables). Similar equation 
can be found in 

The scalar case of ( p4D represents a scalar KP hierarchy in the form of 
superposition formula for three Backlund transformations. To simplify the 
derivation of equations of KP hierarchy from this equation, it is possible to 
take the first order of expansion in terms of the parameter ai, performing, in 
fact, a transition Ai — >■ Ti — > 1 and thus distinguishing the continuous 
variable x := xi (obtaining, in fact, superposition formula for two Backlund 
transformations) . 

It is easy to check that the choice of three shifts Tj : Xi — Xi + [oj] in the 
formula ( |5^ ) in the first nontrivial order of expansion (aiajak) leads to the 
standard form of three waves equation. 

Using cross-differences of equations (^), it is possible to obtain equations 
for the wave functions ('modified equations') 

J2 ^ukTk {{A,T,^) ■ ^-') Ak = 0. (35) 

(ijk) 
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Equations of this type for the DS case and N-wave equations case were derived 
in 



22 , see also the lattice form in 17 



And finally, it is easy to check using (^) that $ satisfies the equation 

(T2Ai$) • (T2A3$)-^ • (T3A2$) ■ (T3Ai$)-i ■ (TiAg^) ■ (TiA2$)-^ = 1. (36) 

This symmetric equation represents a multicomponent singularity manifold 
equation hierarchy. In the scalar case it reads 

(T2Ai$)(T3A2$)(TiA3$) = (T2A3$)(T3Ai<l>)(TiA2$). (37) 

We would like to remind here the main points of the concept of Combescure 
symmetry transformation. The term Combescure symmetry transformation 
originates from geometry, where it is used for the transformation of surface, 
keeping all tangent vectors parallel |^|. For the integrable Darboux system 



describing triply conjugate systems of surfaces this transformation has a very 
simple interpretation in the contexts of theory of integrable systems, and in 
fact a general class of integrable equations ('modified' equations) possesses 
this type of transformation (though the geometrical interpretation in this 
case may be missing). Our concept of the Combescure transformation is 



inspired mainly by integrability, not geometry (see |TT|). 

Integrable equations for the wave functions (i.e. modified equations and 
singularity manifold type equations) possess the symmetry transformations 
corresponding to the freedom of choice of the wave function for the same 
CBA function. In spectral terms this symmetry is described by the change 
of weight functions p(A), ?7(/i). In terms of equations themselves this sym- 
metry is represented using the Combescure invariants. The solutions are 
connected by the Combescure symmetry transformation, if Combescure in- 
variants for them are the same. The Combescure symmetry group consists of 
two subgroups ('left' and 'right'), corresponding to the change of one weight 
function (/(A) or /(/i), see notations (^TJ); in our case /(A) and /(/i) are 
matrix-valued functions, multiplication by /(A) is from the right, by /(/i) is 
from the left, and here we will use the terms 'left' and 'right' Combescure 
invariants according to the rules of matrix multiplication) . Singularity mani- 
fold type equations possess three types of invariants: left, right and complete 
invariants. Modified equations possess only one type of invariants (left or 
right) and are invariant under opposite Combescure subgroup. In this sense 
the basic hierarchy is invariant under the action of all Combescure group. 
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5 Generalized DS hierarchy 

Now we will consider the 2-component case of multicomponent KP hierarchy 
in more detail. In this section we derive the DS hierarchy and the hierarchy 
of corresponding linear problems, the mDS hierarchy and linear problems, 
and in the end we present the compact form of the hierarchy of Combescure 
transformation invariants for the mDS hierarchy. To perform all the deriva- 
tions in terms of hierarchy, we use the formulae containing two continuous 
times (xi := X(i)i and X2 =: X(2)i and discrete shifts [ai] and [02], the first 
acting on X(i), the second on X(2). These shifts change all the times of the 
hierarchy, and the expansion in oi, 02 will generate all the hierarchy in terms 
of partial differential equations. 

To derive DS hierarchy and linear problems, we will use the following 
special cases of the formula ( pSD 



— Aix(A,/i,x) + iri(A)x(A,yU,x) - iri(/i)Tix(A,/i,x) 



= Xi(/W,x)riXi(A,x), 

— A2X(A, /i, x) + K2{X)x{.\ x) - K2{n)T2x{K IJ', x) 



(40) 



= X2(/i,x)T2X2(A,x), 

dix{\ At, x) + (Ki(A) - Ki(/i))x(A, /i, x) = x) xi(A, x; 
(92x(A, /i, x) + (i^2(A) - K2(/i))x(A, /i, x) = X2(/U, x) X2(A, x; 



(41) 

0, (42) 
), (43) 



where we have introduced shift operators and difference operators 



7'ix(A,/i,xi,X2) = x(A,/i,xi + [ai],X2), 
T2x(A, /i, xi, X2) = x(A, fi, xi, X2 + H), 
Aix(A,/i,xi,X2) = (Ti - l)x(A,/i,xi,X2), 
A2X(A,Ai,xi,X2) = (T2 - l)x(A,/i,xi,X2). 
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These operators can also be represented as differential operators 



Ti = exp ^ — <9(i)„ 



n=l ^ 



T2 = exp ^ — <9(2)n • 



\n=l 



The coefficient of expansion of Ti in ai, T2 in 02 are given explicitly in terms 
of Schur's polynomials . . . , %^), p;v(%, • • • , %F) 

The L-operator for the DS hierarchy is obtained by taking (|42|) at A = O2, 
(^31) at A = Oi (and then passing on to the wave functions (pT])) 



82-^1 = w^2, = Xi2(x); (44) 

while the hierarchy of A operators is obtained by taking respectively (^) at 
H = O2, (|41|) at ^ = Oi 

Ai 

—^2 = mTi^i , 

= t;T2^2 • (45) 

02 



One can also obtain dual L and A operators, taking (^OD, ( ^2[ ) at = O2, 
(IT]), (HD at A = Oi. 

It is possible to get DS hierarchy through compatibility conditions, but 
we prefer to extract it directly from the basic formulae (|40|) , (^T]) , (^), (|43|) . 

Taking the difference of equations (^Ol), ( ^2|) at A = Oi, ;U = O2, and also 
using similar combinations, one obtains the system of equations 

(l/ai(ri - 1) - di)u = (Ti - l)xuu, 
(l/ai(l - Tf 1) - d,)v = (T,-' - l)xnv, 

(l/a2(T2 - 1) - d2)v = (T2 - 1)X22V, 
(1/02(1 - T,') - d2)u = (T2-I - 1)X22U, 

dlX22 = uv, 

82X11 = uv. (46) 
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This system of equations represent a compact form of DS hierarchy. The 
expansion of this systems in oi, 02 generates the DS equation and higher 
equations of the hierarchy. It is useful to note that this expansion has the 
structure of recursion relations, expressing the derivatives in higher times 
through derivatives in lower times. 



It is possible to obtain DS hierarchy from the formula (|3^). One could 
start with the three shifts Ti : Xi ^ Xi + [ai], T2 : X2 ^ X2 + [02], T3 : 
Xi ^ Xi + [03] (or T3 : X2 — X2 + [03] for the part of hierarchy with higher 
times connected with the disc 02). Then, taking the first order of expansion 
in Oi, 02 and thus performing a transition from two difference operators to 
derivatives Ai Ti — 1, A2 — T2 ^ 1 and finally writing down 



matrix components of equation (133), one obtains DS hierarchy exactly in the 
form (H). 

The zero order of expansion just gives trivial identities. The first order 
of expansion of the first pair of equations in ai, the second pair in 02 gives 
the DS system 



\ (5(1)2 - ^i) V = -vdixu, 



\ (5(2)2 - u = -ud2X22, 

82X11 = uv, (47) 
where (9(i)2 = a^^' '^(2)2 = d^^- The usual form of the DS system (see e.g. 



[^ ) is a closed subsystem of (|47| ) written in terms of the time t defined as 



dt = 5(1)2 + 5(2)2- 

The second order of expansion, with the use of the DS equations 
gives 

1(5(1)3 — df)u = ud2^di{udiv) — di{ud2^diuv), 

1(5(1)3 - 5i)f = vd2^di{vdiu) - diivd^^diuv), (48) 

1(5(2)3 - dl)v = vd^^d2{vd2u) - d2ivd^^d2uv), 

1(5(2)3 - 52 )m = ud^^d2{ud2v) - d2{udi^d2uv), (49) 
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which is the higher DS system ([^). The reduction u = 1 gives rise to the 
spht Nizhnik-Veselov-Novikov equation for v |25|, [26], while the reduction 
u = V leads to the modified Nizhnik-Veselov-Novikov equation [^. 

To proceed to modified DS hierarchy case, we integrate equations pO|), 
(^TD, (|4^), (^31) over A, /i with some weight functions and obtain equations 
for wave functions 

^$ = ^iTi^i, 
ai 

A2 

^2$ = ^2^2- (50) 

The L operator for modified DS hierarchy can be obtained from the last pair 
of equations, using equations ( ^if ) and dual DS L-operator 



8182^ = (92f/)9i$ + {diV)d2^, (51) 

where f/ = ln\l/i, V = ln\l'2. The hierarchy of A operators arises if one uses 
the second pair of equations (0) to exclude \^2 from the first pair 



A9 T9\l/9 

^$ = 1^292$, W2 = (52) 
02 W2 



The compatibility condition for (pTD, ( p2D gives the modified DS hierarchy 

—Aid2U + {Tid2U){diWi) + {TidiV){d2Wi) + Wi{d2U)AidiV = 
ai 

did2Wi + di{Wid2U), 

—AidiV + Wi{AidiV){diV) - di{WidiV) = 0, 
ai 

Wi = exp(Aif/), (53) 
and for times with the subscript (2) 

—A2diV + iT2diV)id2W2) + iT2d2U)idiW2) + W2idiV)A2d2U = 

0-2 
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d2diW2 + d2{W2diV), 



— A2d2U + W2{A2d2U)id2U) - d2{W2d2U) = 0, 

a2 

W2 = exp{A2V). (54) 
The zero order of expansion of these equations in ai, 02 gives modified DS 



equation ||2^ 



1^2(5(1)2^/ - dfU - {d^Uf) + d,{d2U ■ d,V) = 0, 
9(1)2^ + dlV + {diV) = diV ■ diU, 
ldi{d^2)2V - dlV - {d2Vf) + d2{d2U ■ d^V) = 0, 
di2)2U + dlU + {d2Uf = d2V ■ diU. 



(55) 



(56) 



The hierarchy of Combescure invariants ([Tl|) for the modified DS hierar- 
chy is obtained by representation of the DS solutions u, v through the mDS 
solutions U, V 

u = exp(-Tif/) — exp(V"), 
ai 

Ao 

V = exp{-T2V)—exp{U). (57) 
CI2 



6 r-function for the generalized DS hierarchy 



Now we will analyze the functional equations ([T6| ) for the case N = 2. For 
convenience we present them here explicitly: 



yyi^j Xu{\f^,^+ [a]) 

(6(1) - a(i))xii(A,6(i),x + 
(6(2) - a(2))x2i(A,6(2),x + 

{- — 7^ X22(A,/i,x+ [a 
(6(1) - a(i))xi2(A,6(i),x + 

(6(2) - a(2))X22(A,6(2),X + 



'A 



a(i) 



Xii(A,/i,x+ [b]) + 



VA-6(i)^ 
a])xii(6(i),/i,x+ [b]) + 
a])xi2(6(2),Ai,x+ [b]) = 0, 
f X — a(2) 



(58) 



X22(A,/i,x+ [b]) 



\ A - 6(2) ^ 
a])x2i(6(i),Ai,x+ [b]) + 
a])x2 2(6(2),/i,x+ [b]) = 0, 



(59) 
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Xi2(A,/i,x+ [a]) 



U-^(i), 

(6(1) - a(i))xi2(A,6{i),x + 

(6(2) - a(2))X22(A,6(2),X + 

f^^lx2i(A,/.,x+[a 
(6(1) - a(i))xii(A,6(i),x + 



'A — a 



(2) 



A-6 



(2) 



a])xii(&{i),Ai,x 
a])xi2(6(2),/U,x 
/A-a(i)\ 



Xi2(A,Ai,x+ [b]) + 
+ [b]) + 

+ [b]) = 0, (60) 
X2i(A,/i,x+ [b]) + 



VA-6(i)/ 
[a])x2i(6(i),Ai,x+ [b]) + 
(&(2) - a(2))X2i(A,6(2),x+ [a])x2 2(^2),Ai,x+ [b]) = 0. (61) 

First it is clear that equations ( p8D with a2 = 62 and ( [59| ) with ai = 61 
which coincide with those for scalar KP hierarchy give rise to the following 
expressions for Xii and X22 (see |]T2[) 

1 ri(xi-[A] + [/i],X2) 



Xii(A,/i,x) 



X22(A,yU,x) 



(A-/i) ri(x) 

1 r2(xi,X2 - [A] + [/i]) 



(62) 



(63) 



(A - /i) r2(x) 
where Ti and r2 are the components of the r-function which depend on the 
variables x = (xi,X2). 

Let us consider now equation (|59D at /i = ai and 02 = 62. Substituting 
the expression ( p^ ) into this equation, we obtain 

ri(xi + [61], X2 + [a2])xi2(A,ai,xi + [61], X2 + [02]) = 
ri(xi + [ai],X2 + [a2])xi2(A,ai,Xi + [ai],X2 + [02]). (64) 
For 61 = after the change X2 + [02] ^2, one gets 

ri(xi,X2)xi2(A,ai,Xi,X2) = ri(xi + [ai], X2)xi2(A, 0, Xi + [ai],X2). (65) 
Denoting the r.h.s. of this equation as fi2(A, xi + [ai], X2), we rewrite it as 

fi2(A,Xi + [/i],X2) 



Xl2(A,/i,Xi,X2) 



r^Xi,X2j 

Then the substitution of (|66|) into equation (|6^ ) with 
gives rise to 

fi2(a2,xi + [ai] + [6i],X2+ [02]) 



(66) 

ai and A = 02 



ri(xi + [ai],X2 + [62])- 
r2(xi + [ai],X2 + [62]) 



ri(xi + [ai],X2 + M) 

fi2(62, Xi + [ai] + [61], X2 + [62]) 

ri(xi + [ai],X2 + H) 



(67) 
(68) 
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This equation implies that 



Ti2 (a, xi , X2) = -, 7 (69) 

p(Xi,X2) 



and 



r2(Xi,X2j , , 

pxi,X2 =^ 7- (70) 

ri(xi,X2) 

where T12 is an arbitrary function. Substitution of (|69D into (0) gives 

/\ N ri2(A,xi + [/i],X2 - [A]) . . 

X12 A, /i, xi, X2 = ^ . (71) 

^2(Xl,X2) 

Now we proceed with equation (^). Substitution of the expression (|63| ) into 
this equation with = 02, cti = 61 gives 

X2i(A,a2,Xi + [ai],X2 + [62]) = 

r2(xi + [ai],X2 + [02]) , - 

— — ——^ — 77-^X21 A, 62, xi + ai , X2 + 02 72 

r2(xi + [aiJ,X2 + [O2J) 

After the change xi + [ai] —>■ xi and the choice 62 = one gets 

T"2(x)x2i(A,a2,x) = r2(xi,X2 + [a2])x2i(A, 0, Xi, X2 + [02]). (73) 
From this equation one concludes that 

M N r2i(A,xi,X2 + [/i]) , . 

X21 A,/i,x = — , (74) 

T-2(x) 

where f2i(A,xi,X2) = r2(x)x2i(A, 0, x). 

Substituting now with /i = 02, A = ai, 62 = we obtain 

n(xi + [ai],X2) . , r 1 N n(xi + [6i]^X2) ^ , r, ^ ^ 

— -p- rr2i(ai,xi + ai ,X2) = — —r— 7X21 (61, Xi + pi ,X2). 

T2(Xi + [aiJ,X2) r2(Xi + [6iJ,X2) 

(75) 

Hence 

~ / N T-2i(xi - [ai],X2) , . 

T-2i(ai, xi, X2) = — (76) 

P(x) 

and 

«x)^^. ,77) 
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Substitution of ([76| ) and (^) into ( |7^ gives 



,,.(A,,.x) = ^li^^i^H^iSi±M). (78) 

^1 X 



To establish the relation between ti and T2 we consider again equation ([75|) . 
Due to ( [76D one gets 

ri(xi + [61], X2) ri(xi + [ai],X2) 



r2(xi + [&i],X2) r2(xi + [ai],X2)' 
Since ai and 61 are arbitrary parameters, one concludes that 

n(x) 



(79) 



T-21X 

In the same way, it follows from ( |60D that 



/(X2). (80) 



Thus 



/(xi). (81) 



^ = const. (82) 
r2(x) 

Having in mind the formulae (|62|) and (|63D , one can put ti = T2 without the 
loss of generality. 

So for the 2-component KP hierarchy one has 

V (\ n - 1 r(xi - [A] + [/i],X2) , . 
Xii(^A,/x, xj - — — , (8dJ 

(A - fi) r(x) 

V rA /, ^) - 1 r(xi,X2- [A] + M) 

X22lA,/x, xj - — , (^84J 

(A - ^) r(x) 

,.,(A.,,x)^ ""<"- + Mf -M' . (85) 

r(x) 

r(x) 

At this step r, ri2 and r2i are arbitrary functions. Substitution of the expres- 
sions (|83| - |86|) into the general functional equations (^8| - |6l| ) gives rise to the 
four bilinear functional equations for these functions which represent them- 
selves the addition formulae. We will write down them explicitly for the 
multicomponent case later. 
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7 r-function and addition formulae for n-component 
KP hierarchy 

It will be shown in the next section that the objects of multicomponent KP 
hierarchy depend not only on continuous times, but also on discrete times 
(indices). The addition formulae arising from the discrete dynamics will be 
considered in the context of Toda lattice hierarchy. Here we derive only the 
addition formulae following from the basic formula (|T6|) defining the dynamics 
for continuous variables. 

For the general multicomponent KP hierarchy the system ( ]T^ ) looks 
rather complicated. But its analysis is simplified considerably due to the 
special structure of equations (ITO). It isn't difficult to see that for two fixed 



values of a and (5 , taking = h^^ ^ a, '-^ ^ P, the system reduces to four 
equations of the DS hierarchy (p8|-|6l|) for the variables x„ and x^. Hence, in 
these variables we have the formulae 

^ (\,,.A- ^ r(xi,...,x„ - [A] + [/i],...,XN) 
(A - /i) r(x) 

N Ta/3(Xi, Xq + [/i], X^ — [A], Xn) , //5^/c'7^ 

( A, /i, x) = — ^ (a 7^ /?) • (87) 

r X 



Here indices a and (3 are arbitrary. Therefore the formulae (87) take place 
for all a and (3. 

To compactify these formulae we will introduce the shift operator Tq(A) 
(a = l,...,iV) 

T„(A)/(x) = /(xi,...,x„ + [A],...,xn). (88) 
Using such operators, one can present the formulae ( ^7] ) in the form 

X..(A..,x) = -^^<^iE»i) a = l,...,N. 
A — /i r(xj 

(\ \ 1 T^(/i)T^~^(A)r^^(x) 

Xa/3(A,/i,x) = — a^fJ. (89) 

A — |U r(x) 

In the special cases A = or /x = these formulae reduce to the standard for- 



mulae for the Baker- Akhiezer function ||^, |TJ], |jT6|. Separating the diagonal 
case a = P and offdiagonal one, one gets 

(a, - /i)(A - 6„)r(x + [b])T„(/i)T-i(A)r(x + [a]) + 
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and 



(6, - m)(A - a,)r(x + [a])r„(^)T-i(A)r(x + [b]) + 

(6, - a,)(A - ^)r,(/i)r-i(6„)r(x + [b])r„(6„)T-i(A)r(x + [a]) + 

(/i - A)(/i - 6«)(A - ba) 51(^7 ~ s)^" (/^)^7~^ (^7)^07 (^+ [b]) X 

T,(6,)Tji(A)r,^(x+[a])=0 (90) 



(/i - a,)(A - 6;3)r(x + [b])T,(^)T^-i(A)r,^(x + [a]) + 
ifi - 6„)(A - a;3)r(x + [a])T,(/i)T^-i(A)r,^(x + [b]) + 

- a^){\ - bp)TMTa\KM^ + [h])TMTp\\K/3{^ + [a]) + 
(/i - 6a)(A - 6/3) ^ (b^ - a^)T„(/i)T^^(6^)rQ^(x + [b]) x 



T,(6,)Tji(A)r,^(x+[a]) = 0. 



(91) 



Equations ( pOj ) and ( pij ) represent the extension of the simplest addition 
formula for the r-function in the scalar case (see |jl2[) to the multicomponent 
case. They are related to the generalization of the bilinear relations for the 
Pliicker coordinates to the multicomponent case. On the other hand they 
could be connected with the possible multidimensional extensions of the Fay's 
trisecant formula ETI . 



More general addition formulae for the r-function follow from the general 
functional equation for the function Xai3 which corresponds to the choice 



N 



(^l(x)^2 ^(x))a/3 = 5al3 n 



k=l 



^ — O-ak 
U-bak' 



(92) 



8 r-functional substitution and closed 1-form 

The formulae (|83|), (|8^), (|85|), ([86| ) indicate that it is possible to preserve in 
multicomponent case the general formula connecting the r-function and the 
CBA function in the scalar case. It would be more consistent first to define 
the r-function through the CBA function in terms of the closed one-form, and 
then to prove that this definition implies the general formula. On the other 
hand, one may treat the general formula as the definition and then prove the 
correctness of this definition reducing it to closed 1-form and showing the 



26 



equivalence of 'local' and 'global' formulae (that implies unique existence of 
the r-function, at least in some neighborhood of g = 1) . In this paper we 
would rather follow the second way. 

So we define the r-function by the relation 

x(A,/i,^7)=xo(A,/i)%^. (93) 

-rig) 

Here Xo(A, fi) is just the function (A — /x)^^ extended as 1 for A, /x in different 
Di, and the function g(x,fi) is an extension of the function to the case of 
disconnected G 

Xo A,/i = V f ^^ n ^ n ■ ^'^ 

[ 1 ioT X e Di, fi e Dj,i ^ J . 

The function g(x^^) ( in fact, it is a set of loop group elements gii^i') de- 
pending on the parameters A, fj,) is the ratio of two terms 

g{\) 

g{x,t.) = — , 
gill) 

r \ \ V for A G A /n.x 

3W^i^) = \ 1 for A ^ A. ^^^^ 



To extract the formulae (|B|), (||), (|H|), (|§ from the formula (|53|), one 
should recall that the loop group defining the multicomponent KP hierarchy 
can be decomposed in the following way 

N 

giii^) = gto,M X z/"% Y.ni = 0, (96) 

1=1 

where the first term defines the dependence on continuous times and the 
second - on the discrete times (indices). If A, /x belong to the same Di, using 
explicit parametrization of the loop group elements, we obtain exactly the 
formulae (|5BD, (Pll). For A G -Dfe, /i G Dp, k ^ p, g{\,^) has a factorization 

g{\,ii)k[v) = X J^, 

u 

I \ ^ -1 

z/ — /i 

fi'(A,M)»('^) = 1, i^k,i^p. (97) 
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Using these formulae, we obtain the formulae (^), (p6D, with the compo- 
nents of r-function interpreted as a r-function with shifted index variables 
(general index zero is preserved). For example, in the two-component case 
the function T12 corresponds to r(ni — 1, n2 + l). To express the CBA function 
with zero indices, we need to know the set of functions r, Tij, i ^ j. Then for 
other admissible index vectors the CBA function and the r-function can be 
found explicitly through the determinant formulae |0. Thus, the discrete 
dynamics is exactly solvable in terms of initial data (CBA function with zero 
indices or the set of r-functions r, t^, i ^ j). 

Returning to our original purpose, we will extract the definition of the 
r-function in terms of the closed 1-form from (^ and then prove that both 
definitions are equivalent. Our derivation will mostly repeat the derivation 



in the scalar case given in |]T2|, but we would like to remind it due to the 
importance for the following consideration. 

Differentiating (p3|) with respect to A, one obtains 



1 r 6t (g X 5'(A,/i)) 1 d 

c/z/ = — (A -/i)xaa(A,/i;5f) , (98) 



T{g) I S\nga{iy) v - \ dX 
or for A = /X 

1 / Srig) 1 

TT r Ti rdi' = Xraa{\X;9), (99) 

T{g) J din ga{iy) u - X 

where Xra/3(A, /i; g) = XaisiX, fj,; g) — 5ap{X — is for 'regularized'. This 

equation can be rewritten in the form 

brig) 1 / 1 

-du = <p XTaa{^,J^;g) rdiy. (100) 



T{g) I 6\ngaii^) u - X I ' ' u - X 

Relation ( p.00| ) implies that the linear functionals ~:;j^ ^^sing^ Xri^, ^'1 9) 
are identical for arbitrary variations y analytic outside each Di and decreas- 
ing at infinity. Thus for y belonging to this class of functions the variation 
of the r-function is given by 

- 6\nT{g) = f Xraa{j^,J^;g) . du . (101) 



a=l 
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This expression defines a variational 1-form defining the r-function. It is easy 
to prove using the identity ([l|) that this form is closed . Indeed, according to 
(|I]) (here we are using scalar notations, it is not difficult to extract equations 
in terms of components Xa/si^, f^', g)) 

^x{>^, /w; ^) = / x{^, f^; 9)^^7-lx{\ 9)dy , 
JdG g[v) 

^Xr(A, A;5() = / Xr{.^A]g)—r^Xr{.\j^]g)di'. (102) 

JdG g{iy) 
So the variation of the ( |101| ) gives 

-6'lnT{g)=ff xMX-g)x,{X,u-g)^-^^^dudX. (103) 

The symmetry of the kernel of second variation with respect to A, v implies 
that the form ( |101[ ) is closed. 

So the formula ( |101| ) gives the definition of the r-function in terms of the 
closed 1-form. For the standard KP coordinates 

ga{X) ~h A" 

This formula allows us to obtain a closed 1-form in terms of dx^n 



N oo Qn 

rflnr(k,x) = XI ^Xraa(i^, '^,k,x) 



dX(a)n ■ (104) 

u=0 



To prove that the 'local' definition in terms of the closed 1-form implies 
the global definition (^), it is sufficient to show that logarithmic derivatives 
of the l.h.s. and the r.h.s. of ( p3| ) with respect to A, /i and variational 
logarithmic derivatives with respect to Sga/ga are equal in virtue of ( |101j ) 
and generalized Hirota bilinear identity (|I]). Then the l.h.s. and the r.h.s. of 
( |93[) could differ only by the factor, and the normalization of the function x 
implies that this factor is equal to 1. 

Here we will prove only the equality of logarithmic variations, the proof 
for derivatives is analogous (see []T2|). 



29 



The logarithmic variation of (^3]) gives 



6 (lnx(A,/i;5()(A - /i) - \nT{g x g(^x,,,)) + ^^r{g)^ 




{Xt{j^, ^-,9^ 9(x,^)) - Xt{i^, 9)) 




(105) 



(we have used ( |102| ) and ( [101| )). After the use of the determinant formula 



to transform the second term, one concludes that the variation is equal to 
zero. 

9 Generalized 2D Toda lattice hierarchy. Ba- 
sic equations. 

We will consider here for simplicity only scalar 2D Toda lattice hierarchy. 
It is generated by the generalized Hirota identity with the domain G which 
consists of two disconnected regions, i.e. in some sense it is a two-point KP 
hierarchy. But it also possess some fundamental special features and can be 
treated as a primary object with respect to scalar KP hierarchy. 
Without loss of generality one can choose G as 



where Dq and D^o are unit discs around the origin and infinity respectively. 
The two-component KP hierarchy is generated by the pair of functions (7in, 
(7out defined on the boundary of Dq, Doo respectively, having zero general 
index (let their individual indices be n, —n), gout is analytic and have no 
zeroes outside the unit circle, ^^in inside the unit circle. It is easy to see, that 



X(A',/^') x(A',/i) 
x(A,/i') x(A,/i) 



(106) 



x(A',/i') 



G = DoUD, 
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dynamics defined by the pair of these functions coincides with the dynamics 
defined by one function for both boundaries 



g = A"(7in^7out (107) 



oo 



up to some trivial dynamics, defined by a function analytic m G = DqU D 
{g = gin, \e Do; g = gant, A e D^). 

As we know, arbitrary function on the circle having no zeroes (arbitrary 
loop) can be represented in the form ( 107 ). So in the case of Toda lattice 



hierarchy the objects of the theory nontrivially depend on the whole loop 
group (i.e. on discrete time, positive and negative continuous times), r- 
function in this case is a functional of arbitrary g on the circle. Discrete 
dynamics is correctly defined, because in this case A" is a function having in 
G equal number of zeroes and poles. 

The explicit parametrization of the loop group element g in terms of times 
is given by 

5((n,x, z/) = 5((n,x+,x_,i/) = z/^exp ^i+)k'^'^ + ^{-)k'^''^ ■ (108) 

The compact form of 2D Toda lattice hierarchy in terms of the function 
x(A,/i, /c,x) can be obtained from (|l]) using special combinations of group 
elements gig2^ with simple analytical properties. To define dynamics with 
respect to continuous times, let us take 

^o(-):=^7i(%2 ^ = ^^7357' 
a+,b+ e Do, a_,b^ G D^, 

go{u) = g+ X ig_)-\ g+ = -±, g_ = . 109 

z/ — 0+ V — a_ 

Substituting the expression for gig2^ into the generalized Hirota bilinear 
identity (|lD, one obtains 

^o(/^)x(A, ^^,gx g+) - ^o(A)x(A, /x, ^ x g_) + 
Res(5(o)|6+ x(A,&+,5' x g+)xib+, fi, g x g_) + 
Res(^o)|b_ x{Kb_,g x g+)x{b^, fi, g x g_) = 0. (110) 
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The equations of Toda lattice hierarchy are generated by equation ( |110| ) with 
6+ = 0, b- = oo, or, in terms of go 

goiiy) = 9i{^)92» = (1 - -)(i - -) = exp ff: -{{-y + {-y)] m) 

|a_i_| < 1, |a_| > 1 . 

Substituting into equation ( |110D parametrization of loop group elements in 
terms of times, we obtain 

(1 - — )(1 - — )x(A,/i,n,x+ + [a+],x_) - 

(1 - — )(1 - /i, n, x+, x_ - [a_]) + 

-a+ x(A, 0, n, x+ + [a+], x_)x(0, /i, n, x+, x_ - [a_]) 

X{\ oo, n, x+ + [a+], x_)x(oo, /i, n, x+, x_ - [a_]) = 0, (112) 

where x(A, cxd) = ;\;(A, /i) x /i|^=oo- In terms of CBA and BA functions ^'(A, /i) 
and Vo(A), ^o{li) ^oo(A), ^/'oo(/i) we have 

/i, n, x+ + [a+], x_) - /i, ra, x+, x_ - [a_]) = 
a+?/'o(A,n,x+ + [a+], x_)^o(/i, ^, x+, x_ - [a_]) + 

— ?/^oo(A,ri,x+ + [a+],x_)?/'oo(yU,?^,x+,x_ - [a_]). (113) 

(X — 

To get equations defining the discrete dynamics, we take 

9M=9M92\y) = \ (114) 

or 

9^{y)=9M92\^) = ^- (115) 
In the first case the equation for the function x reads 

-X(A, /i, ri + 1, x+, x_) - \x{\ yu, x+, x_) = 

X(A, 0, n + 1, x+, x_)x(0, /i, n, x+, x_), (116) 
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and in the second case 

/iX(A, /i, n, x+, x_) - Ax(A, /i, n + 1, x+, x_) = 

X(A, oo, n, x+, x_)x(oo, /i, n + 1, x+, x_). (117) 

In terms of CBA and BA functions one gets 

^p{X, /X, n + 1, x+, x_) - 'ip{X, /i, n, x+, x_) = 

^o{\, n + 1, x+, x„)^/io(/i, n, x+, x_) 

^p{X, /i, n, x+, x_) - ^Ij{X, /i, + 1, x+, x_) = 

V'oolA, n, x+, n + 1, x+, x_). (118) 

Equations ( [L13| ), ( |118| ) define the generahzed 2D Toda lattice hierarchy. The 
discrete dynamics can be resolved in both directions explicitly (i.e. x(n + 1), 
x{iT- — 1) can be expressed through x{^)) using the determinant formula 



10 Toda lattice hierarchy: 2DTL, 2D Volterra 
chain, 2DTL singularity manifold equa- 
tion 

At this point we would like to derive the simplest equations of generalized 2D 
Toda lattice hierarchy. They contain three variables: one discrete variable n 
and two continuous variables (we will denote them x+ and a;_). 

We will derive three different equation corresponding to different levels of the 
tower of equations of generalized hierarchy: the main equation (the standard 
two-dimensional Toda lattice equation), the modified equation (i.e. equation 
for Baker- Akhiezer function) and singularity manifold type equation (arising 
here as an equation for the CBA function). It is interesting to note that the 
second and the third equation (not to mention the first) are already known 
in literature in different contexts (see |]TP|], [0, PU[). 



The basic equations for our derivation are equations (|116|) , ( |1171) and the 
zero order terms of expansion of ( p.l2| ) in a+, l/a_, i.e. 

(^ + ^--)x(A,/.,n,x^,x_) = 
X(A, 0, n, x+, x_)x(0, /i, n, x+, x_), 
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d 



"dx. 



+ A - /i)x(A,/i,n,x+,x_) 



X(A, oo, n, x+, x_)x(oo, /i, n, x+, x_). 
In terms of CBA and BA functions instead of ( |119| ) we have 
d 



(119) 



dx^ 
d_ 
dx- 



^(A, /X, n, x+, x_) = V^o(A, ra, x+, x_)?/^o(/x, x+, x_), 
'0(A, /i, n, x+, x_) = ?/'oo(A, n, x+, x_)Vioo(/i, n, x+, x_) 



(120) 



Combining ( |116| ), ( |117| ) and ( p.l9| ), one obtains 



8 1 1 

X(A, 0, n, x+, x_)(x(0, /i, n, x+, x_) - x(0, /i, n - 1, x+, x_)), 

(^^ + A)x(A, /i, n, x+, x_) - Ax(A, fi,n + l, x+, x_) = 

X(A, oo, n, x+, x_)(x(oo, fi, n, x+, x_) - x(oo, /i, n + 1, x+, x_))(121) 

Taking the first equation at A = oo, /i = 0, the second at A = 0, yU = oo, one 
obtains the following equations 



d+u{n) = {U{n) — U{n — l))u{n), 
d_v{n) = {V{n) - V{n + l))v{n). 



(122) 



Here we have introduced the notations 

u{n) = x(oo,0,n,x+,x_), 
v{n) = x(0,oo,n,x+,x_), 
Uin) = Xr(0,0,n,x+,x_), 
Vin) = Xr(oo, oo, n, x+, x_), 

and we omit arguments x+,x_ common for all functions. 

Equations ( [L19| ) taken respectively atA = oo,/i = oo; A = 0,/i 
rise to the following relations 

d+V{n) = u{n)v{n) — 1, 
d-U{n) = u{n)v{n) — 1, 



give 



(123) 
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which together with ( |122| ) form a closed system of equations. 
There is also an additional relation implied by (|116|) , ( |117| ): 



u{n + l)v{n) = 1. (124) 

This formula shows that the functions u{n) and v{n) are not independent. 

The system of equations (|122|) , (|123|) , (|124|) represent some special form of 
classical 2D Toda lattice equations. To show it, let us derive more standard 
forms of Toda lattice from this system. First, it is easy to show that 

c}+9_ \nu{n) = u{n)v{n) — u{n — l)v{n — 1), 

(9+9- \iiv{n) = u{n)v{n) — u{n + l)v{n + 1). (125) 

Thus for the function h{n) = u{n)v{n) we have an equation 

d+d^\nh{n) = 2h{n) - h{n - 1) - h{n + l), (126) 

which is 2D Toda lattice equation in the Darboux form. 

Since, according to (|124D, the functions u{n) and v{n) are not indepen- 



dent, we can write an equation for one function 

r^ -, / N uin) uifl — l) / „n 

d+d^ Inw n = , \. - , , • 127 
^ ' u{n + l) u{n) ^ ' 

In terms of Qin) = In 'u(n) we have 

d+d.e{n) = exp{e{n) - e{n + 1)) - exp(0(n - 1) - e{n)). (128) 

This is just a standard form of the 2D Toda lattice equation. An original 
bilinear form, introduced by Toda, follows from addition formulae for the 
r-function, which will be given in the next section. 

Now we proceed to the modified Toda lattice equation. In this case we 
will use the equations for generic wave functions, which are given by ( |11^ ) 
and ( |120| ) integrated over A and /i with some weight functions p(A) and ?7(/i) 
(in fact, each of them in matrix notations is a pair of functions defined on 



'in' unit circle and 'out' unit circle, see also (|3l|) ) 



$(n) - $(?2 + 1) = ^o('^ + l)^o(^), 

<l>(n)-$(n-l) = ^oo(ri-l)^oo(n) (129) 
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(we have omitted common arguments x+,x_) 

Combining these equations, one obtains a hnear problem for the modified 
Toda lattice equation 



= W+{1 - T)$, 



= 1^_(1 - T-i)<l>, W_ = (130) 

The compatibihty condition gives the modified Toda lattice equations 

d_W+ = W+{1-T)W_, 

d+W. = iy_(l - T-^)W+. (131) 



These equations quite recently appeared in literature |jT9|, [T^ under the 
name of two-dimensional Volterra chain. It is interesting to note that in [|18 
they arose in the geometrical context. 

Proceeding further, one can obtain an equation for the function $, using 
linear operators of modified Toda lattice to express the functions 1^+, 



via the function $ and then substituting the final result to equations ( |131| ). 
Both equations give rise to the same equation 

d+d-^ + 9+$ ■ ( ^^ ^ + -—^ — -] = 0. (132) 

This equation was also introduced in different contexts in [|T^, [|T^]. Another 
interesting feature of equation ( |132| ) could be its connection with the Painleve 



analysis of the 2DTL. Making comparison with one- dimensional Toda lattice 
pO[| , and recalling analogous equation for the scalar KP case, it is natural 



case 



to suggest that the equation ( |132[ ) represent a singularity manifold equation 
for the 2DTL. 



Equation (|132|) generates 2DTL and two-dimensional Volterra chain in 
different ways: first, in terms of Combescure invariants, and second, in terms 
of expansion in A, fi of the function ip{X, /i, n, x+, x_) which satisfies equation 

The Combescure invariants for ( |132| ) are given by the formulae 

W. = (133) 



36 



(we have written down only 'left' invariants). For two-dimensional Volterra 
chain the Combescure invariants look like 



u = W. exp{-A-\W. + TW+)), 

V = W+exp{A'\W.+TW+)). (134) 

Substituting ([133D to (|13|), one obtains expressions of the full Combescure 
group invariants for singularity manifold equation ( |132D . 

At last, we present the formulae for the hierarchies of two-dimensional 
Toda, Volterra and 2D Toda singularity manifold equations. The derivation 
is very similar to the case of DS hierarchy, so we will omit it. 

The 2DTL hierarchy can be represented as an equation 

T^T^'u + ^Tu = T^uT-\ f 1 + \ ] , (135) 

where a convention for shift and difference operators is standard 

r+ : x+ ^ x+ + [a+], 
T_ : x_ ^ x+ + [aZ^]-, 

A+ = r+ - 1, 

A_ = (T_ - 1). (136) 

The 2D Volterra hierarchy reads 

a^A_W+ + {T_W+){TW^) - {T+WJ)W+ = 0, 

l/a+A+W^_ + {T+W-){T-^W+) - {T_W+)W- = 0, (137) 



where Wm = + ° W - 

Finally, 2DTL singularity manifold equation hierarchy looks like 

{r_(T+ - 1)$}{T(T_ - r"^)$}{T+(l - T"^)$} = 

{(T+ - l)$}{T+(r_ - T"^)$}{T_(r - 1)$}. (138) 

This equation possesses a symmetry ^ T_, T_ ^ T ^ T^^. Though 
this symmetry is not obvious but it can be easily verified. 
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11 Addition formulae for Toda lattice hierar- 
chy 

The r-functional substitution 

X(K,.9)- ' "'y' , (139) 
X- 12 T{g) 

X-u 

fi'(A,/.) = , 

fl — U 

to ( P-lOp , ( p.l8| ) leads to addition formulae for the r-function, the expansion 
of which gives the equations of the Toda lattice hierarchy in terms of the 
T-f unction. 

The simplest addition formulae follow from ( |118|) 



-^{g X 9dg{x,^))r{g) - — r-^ — -T{g x g(x,,,))r{g x g^) + 



j^r{g X gdgixfl))r{g x ^(o,;,)) = 0. (140) 

This is a functional form of addition formulae, containing four compo- 
nents characterized by the position of A, fi with respect to the unit circle. 
Using parametrization of loop group elements in terms of times, we obtain a 
standard form of addition formulae. It is convenient to write down first the 
factorization formulae for the function g(\,fj,), representing it as a product of 
functions analytic inside and outside the unit circle for different positions of 
A and /x. For A, /x both belonging to Dq (Z^oo) the factorization is trivial, 
because g{\^^) is analytic outside (inside) the unit circle. For A, /i in different 
domains it is given by 

g{X+,fi^) = 5'd5'(A+,0)5'(oo,/i_), 
g(X-,fi+) = fi'd ^fi'(A_,oo)5'(0,/i+), 

gd{i^) = (141) 

here the subscript +, - just indicates the position inside or outside the unit 
circle. 

In terms of time variables, the action of loop group elements represented 
in the formula (^30) (taking into account the factorization formulae) looks 
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like 

9ix+,^l+) ^ x+ + - [A+], 

9{x+,t,-) + l,x+ - [A+],x_ + [/il^], 

9{x-,ii+) ^ n - l,x+ + [/i+],x_ - [Al^]. (142) 

Substituting these formulae to ( |140| ) , we obtain a standard form of addition 
formulae for the Toda lattice hierarchy. It occurs that off-diagonal compo- 
nents of ( |14CI| ) give the same formula with interchanged A, /i, so we have three 
addition formulae 

.{x+ + [,,],x.) {x+,x. + [X-^) _ ix+ + [,,],x-) ix+,x^ + [X-^]) _ 

(A - ^)t-^-+.--)^^-++M,--+[a-^])^ 

|A| > 1, < 1, (143) 

\ + {x+ + [X],x-) _ (x+ + [tj.],x-) {x+ + [X],x-) _ 

'n f^'n 'n+1 — 

(A-/^)ri:r-M^^^f^^-'''^'^-\ 

|A| < 1, < 1, (144) 

^(x+,x. + [^,-^) (x+,x^ + [X-^]) ix+,x_ + [,i~^]) (x+,x_ + [X-^]) _ 

^'n 'n-1 P'n-l ~ 

(A-/i)er-M^-^-^'^"^+''"^\ 

|A| > 1, > 1. (145) 



In fact, we have used only the first equation of the pair of equations ( |118| ), 
but the second equation gives exactly the same set of addition formulae. 

The formula (|143|) is a well-known simplest addition formula for the Toda 
lattice hierarchy. The formulae (|144|) , (|145|) do not shift negative (positive) 
times and it is connected with the Backlund transformations for the KP 
equation. 

To get more generic addition formulae, one could use the formula ( p,lCI| ). 
To get more symmetric expressions, we will use a special version of this 
formula corresponding to the choice 

9o{iy)-=9i{'^)92\^) = ^^, (146) 
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where we do not specify from the beginning the position of the points Ai, /ii. 
This choice leads to the following functional form of the addition formula 



(/i - r/)(A - Or(^^(^,^)^(A,^))r(^) -{\-ri){fi- 0^{gg{x,^,))T{gg^^.r,)) + 
(A Origgiii,,Mgg^x,,)) = 0. (147) 

The formula ( |147] ) contains 16 components depending on the position of four 
parameters with respect to unit circle. It is convenient to introduce the 
function showing the position of each parameter 







< 1 




A 


> 1 



^A= \ (148) 



Then the addition formulae take the form 

(M-.)(A-Ora»"'etlT™"- 

(A - - Oet!L1;"^"'et^*"«» = 0, (149) 

where 

_ f (x+ + [A],x_) |A| < 1 

i (x+,x_ + [A-i]) |A|>1 

Due to its symmetric form, the formula ( p.49| ) contains only five different 
addition formulae corresponding to all four parameters inside the unit circle, 
one parameter outside, two parameters outside, three parameters outside, 
all four parameters outside. The first and the last case just give the KP 
addition formulae, so we will consider only remaining three cases. First we 
write down the addition formulae for the case when two variables (say, A and 
/i) are inside the unit circle: 

-7]){X- ^)^{x+ + M,X- + [,,-i])^(x+ + [A],x_ + [rl]) _ 
(A - r/)(/i - ^)^{x+ + M,x- + [r^])^(x+ + [A],x_ + [r,-i]) ^ 

(A -f^)iv- = 0. (150) 



The structure of this addition formula is similar to the formula ( |143| ), and in 
fact ( |150| ) is a generalized version of ( p.43| ). 
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When one parameter (say, A) is inside the unit circle, we obtain the 
following addition formula 

(A -fi)iv- e)r(^++W'''-+['^"'])ri^t""^''"'^'^"'^ = 0. (151) 
And the last case (A outside the unit circle) gives 

- r/)(A - ^)^(x++W+M,x-)^x++K],x-+[A-]) _ 

(A ^)r(-++M--+[A-^])r(^/M+K]'--) = o. (152) 

The formulae (|l5l]), ( |152D generalize the formulae (|14|), (|T45|). 
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